Heuristic theory for many-faced 
d- dimensional Poisson-Voronoi cells 



H. J. Hilhorst 

Laboratoire de Physique Theorique, Batiment 210 
Universite Paris-Sud and CNRS 
91405 Orsay Cedex, France 

June 24, 2009 



Abstract 

We consider the (i-dimensional Poisson-Voronoi tessellation and inves- 
tigate the applicability of heuristic methods developed recently for two 
dimensions. Let p n (d) be the probability that a cell have n neighbors 
(be 'n-faced') and m n {d) the average facedness of a cell adjacent to 
an rz-faced cell. We obtain the leading order terms of the asymptotic 
large-n expansions for p n (d) and m n (3). It appears that, just as in 
dimension two, the Poisson-Voronoi tessellation violates Aboav's 'lin- 
ear law' also in dimension three. A confrontation of this statement 
with existing Monte Carlo work remains inconclusive. However, sim- 
ulations upgraded to the level of present-day computer capacity will 
in principle be able to confirm (or invalidate) our theory. 

Keywords: Poisson-Voronoi cell, number of neighbors, two-cell cor- 
relation, Aboav's law 
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1 Introduction 



1.1 General 

Cellular structures occur in a wide variety of natural systems. The examples 
most quoted, but by no means the only ones, are soap froths and biological 
tissues. Cellular systems also serve as a tool of analysis in a diversity of 
problems throughout the sciences and beyond. Many references may be found 
in Okabe et al. [1], Rivier [2], and Hilhorst [3]. 

A popular model of a cellular structure is the Voronoi tessellation. It is 
obtained by performing the Voronoi construction on a set of point-like 'seeds' 
in (i-dimensional space. This construction consists of partitioning space into 
cells in such a way that each point of space is in the cell of the seed to which 
it is closest. A ^-dimensional Voronoi cell is convex and is bounded by planar 
(d— l)-dimensional faces. Two cells are called 'neighbors' or 'adjacent' when 
they have a face in common. The Voronoi construction may thus serve to 
define neighbor relations on an arbitrary set of given seeds. 

In the special case that the seed positions are drawn randomly from a 
uniform distribution, one speaks of the Poisson- Voronoi tessellation. It con- 
stitutes one of the simplest and best studied models of a cellular structure. 
The analytic study of the statistical properties of the Poisson- Voronoi tessel- 
lation in dimensions d = 2, 3 was initiated by Meijering [I] in 1953. Monte 
Carlo results were obtained by several workers in the past decades (see pQ 
for references). The analytical and numerical results of greatest interest in 
dimensions d = 2 and d = 3 have been listed in Ref. pQ. In arbitrary di- 
mension d, a large collection of statistical properties of the Poisson- Voronoi 
tessellations were derived in Refs. [3, [5] . 

One of the most characteristic properties of a cell is its number n of faces, 
and it so happens that this is not a property easy to study analytically. There 
must be at least d+ 1 faces, but there may be any number of them. We denote 
by p n {d) the probability that a randomly picked cell (all cells with the same 
probability) be n-faced. The facedness probability p n (d) has been the center 
of interest of experimental and Monte Carlo work (see [HO [3] for references), 
especially in d = 2 and d = 3. Analytical results for this quantity, however, 
are very scarce. Only the one-dimensional case with p n (l) = is trivial. 
Known results in higher dimension include the average facedness (n)d = 
^2 n np n (d) in dimensions d = 2,3,4. Its values are (n) 2 = 6, (n) 3 = 2 + 
48tt 2 /35 = 15.535..., and (n) 4 = 340/9 = 37.77... [7J. It has been determined 
numerically that the peaks of the distributions are at n = 6 for d = 2 and at 
n = 15 for d — 3; as n increases beyond the peak value, p n (d) decreases very 
rapidly to zero. 
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1.2 Recent work 



Although in general dimension d one readily writes down an n<i-dimensional 
integral for p n (d), the variables of integration are coupled in such a way that 
for all d > 1 this is a true many-particle problem. As a consequence, it 
has not been possible to determine the neighbor number probability p n (d) 
analytically. In dimension d = 2 progress was made, nevertheless, in Refs. [H 
[3], where it was shown that in the limit n — > oo the sidedness probability 
p n (2) has the exact asymptotic behavior 

Pn(2) = ^^[1 + 0(1)], r^oo, (1.1) 

with C = 0.344 347.... The collection of asymptotic results that include (11.11) 
required considerable calculational effort. Subsequently, however, heuristic 
arguments were developed [9] by which at least part of the same results can 
be derived more easily. We extend in this work these heuristic methods to 
higher dimensions. 

The correlation between the neighbor numbers of adjacent cells is usually 
expressed in terms of the quantity m n (d), defined as the average facedness 
of a cell that is itself adjacent to an n-faced cell. Aboav's celebrated 'linear 
law' states that nm n (d) = an + b [10J. It holds within error bars, and with 
system-specific a and b, for many experimental cellular systems. In Ref. [TT] 
it was demonstrated, however, that for the two-dimensional Poisson-Voronoi 
tessellation Aboav's law is in fact a linear approximation limited to small n 
values. The true asymptotic behavior appeared to be 

m n (2) = 4 + 3(vr/n) 1/2 + . . . , n -> oo, (1.2) 

and shows that nm n (2), instead of being linear, has a small downward curva- 
ture. In fact, this deviation from the linear law was known from Monte Carlo 
simulations [T2], [131 [H] . In this paper we investigate how (II. 2p is modified in 
dimension d = 3. 

1.3 This work 

This work continues a series of articles [H |3l [TTJ, [151 HS1 E] that deal with the 
properties of Voronoi tessellations. We consider here the question of what, 
if anything, may be learned in higher dimensions from the two-dimensional 
case. The answer is that the exact methods, as in so many other domains of 
physics, remain limited to two dimensions of space. However, we can apply 
to higher dimensions the heuristic methods that we developed [9] on the basis 
of the exact ones. 

In this way we obtain first, in section [21 we obtain a formula for the 
large-n behavior of p n (d) in arbitrary space dimension d. The argument is 
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a fairly straightforward extension from the two-dimensional case [9]. Then, 
in section [3j we obtain a two-term asymptotic large-n expansion of m n in 
dimension d = 3. Our key results are represented by equations (12.201) and 
(13.21) . that are the higher-dimensional analogs to (II. ip and (11.21) . respectively. 
One of our findings is that Aboav's law fails also for the t/iree-dimensionsal 
Poisson-Voronoi tessellation. In section H] we compare our theoretical result 
for m n (3) to existing Monte Carlo work. We conclude in section [5j 

2 Asymptotic large-n expression for p n (d) 

We consider a <i-dimensional Poisson-Voronoi tessellation with seed density 
p. This density may be scaled to unity but we will keep it as a check on 
dimensional coherence. We select an arbitrary seed, let its position be the 
origin, and are interested in its Voronoi cell, to be referred to as the 'central 
cell'. The central cell has the same statistical properties as any other cell. 
We set ourselves as a first purpose to determine the facedness probability 
p n (d) of this cell in the limit of large n. 

2.1 The shell of first neighbor seeds 

In the two-dimensional case the n-sided cell is known for large n to tend with 
probability 1 to a circular shape. It is natural to assume that in d dimensions 
the n-faced cell will similarly tend to a hypersphere when n —>■ oo. Let R 
denote the n dependent radius of this sphere. Then the central seed has its 
n first-neighbor seeds located, for large n, within a narrow spherical shell of 
radius 2R. We denote the width of this shell by w; approach to sphericity 
means that w/R — > as n — > oo. 

This leads us to the 'shell' model for the first-order neighbor seeds repre- 
sented in figure [H We consider two d-dimensional hyperspheres of radii 2R 
and 2R — w, both centered around the central seed. Let all seeds other than 
the central one be distributed independently and uniformly in space with 
density p. Let now p n (d;R,w) denote the probability of the event - which 
defines the 'shell' model - that (i) the inner hypersphere contains no other 
seeds than the central one; and (ii) the shell contains exactly n seeds. Our 
procedure will be to write p n {d] R, w) as an explicit function of n, d, and the 
two parameters R and w (this is easy). We will then, by means of a heuristic 
argument, express w in terms of R and finally maximize with respect to R. 
The result will be the desired expression for p n (d). 

The probability p n (d; R, w) follows from an elementary calculation and is 
equal to 

p n (d;R,w)=(^e-» v >, (2.1) 
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Figure 1: Schematic representation of the 'shell' model of the first-neighbor seeds. 
Two hyperspheres (of radii 2R and 2R — w) are centered about the origin O of 
(i-dimensional space. A 'central' seed is located in O. Other seeds are distributed 
randomly and uniformly with density p. The picture shows the rare event of having 
(i) n seeds inside the shell; (ii) no seed inside the inner hypersphere. The circle 
of radius R is the approximate boundary of the central Voronoi cell. The asterisk 
denotes the location of contact between the central cell and the cell of seed F2. 
Further explanation is given in the text. 
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Figure 2: Detail of figure Q] in which seed F 2 has moved a little outward and is 
now screened by its neighbors F\ and F3 from being in contact with the central 
cell. 

where V\(R,w) and V^-R) are the volumes of the shell and of the outer 
sphere, respectively. We will denote the volume v d and the surface area s d of 
the o?-dimensional hypersphere of unit radius by 

v d = J/(l)\, Sd = 2^/(1-1)1 (2.2) 

Therefore 

Vx(R,w) = v d {(2R) d -(2R-w) d }, 

V 2 (R) = v d (2R) d . (2.3) 

It will be convenient to work with \ogp n . Using (12.31) in (12. ip we find 

\og Pn (d; R,w) = - log n\ -log [pv d {{2R) d - {2R - w) d }] -pv d (2R) d , (2.4) 

which is exact within the shell model. This expression for p n (d; R, w) will be 
at the basis of what is to follow. 

2.2 Relation between w and R 

We now look for a relation between w and R. Figure [T] is a schematic two- 
dimensional representation of a <i-dimensional situation. It shows part of the 
Voronoi cell boundaries of three first-order neighbor seeds F\, F%, and F3. 
Each line segment in the figure belongs to a (d — l)-dimensional face that 
perpendicularly bisects the line segment joining two neighboring seeds. An 
asterisk marks the short line segment common to the central cell and the cell 
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Figure 3: Marginal situation in between those of figures [T] and EJ The central 
cell has a point contact with the cell of seed F2 in the four-vertex denoted by the 
asterisk. 
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of F 2 . As shown in figure [21 this face is so small that it disappears when 
F 2 , roughly speaking, crosses to the outside of the outer hypersphere: F 2 is 
then be 'screened' by its neighboring seeds F% and F3, When w is small with 
respect to the typical interseed distance I in the shell, then such screening 
will be negligible for a random arrangement of the n seeds in the shell. We 
now relate shell width w to the radius R by imposing that w be small enough 
for this to be the case, but otherwise as large as possible. Figure [3] depicts 
a marginal situation in which the screening of F 2 by Fi and F 3 sets in. 
The idealization consisting in placing Fi and F3 at equal distances from the 
origin and symmetrically with respect to F 2 is good enough for our ourpose. 
Elementary geometry then shows that w, £, and R are related by 



w l , 
- ~ — . (2.5) 



I 2R 

where the symbol ~ denotes proportionality in the limit of large n. We now 
wish to eliminate t from this relation. 

The shell may be considered locally as a flat id — l)-dimensional space if 

(2.6) 

We will asume here, and be able to verify afterwards, that in the limit of large 
n these conditions are satisfied. The typical interseed distance I between the 
first-neighbor seeds in the shell is then easily found. For w <C I the n seeds 
may be considered as distributed on a (d — l)-dimensional hypersurface of 
area 

S = s d {2R) d -\ (2.7) 

They therefore have a (d — 1) -dimensional surface density a = n/S, whence 
it follows that their typical distance i may be defined by 

v d _ x i d - x = Sn-\ (2.8) 

From and §ZB) it follows that 

£ = 2r( -^O . (2.9) 
\nvd-xj 

Comparing finally (12.91) and (12.51) suggests that in the shell model we should 

set 

w = 2R(c d n)~^, (2.10) 

where q is a (not exactly known) numerical constant of order unity. Whereas 
have argued above for the validity of (12 . 10[) in the limit of asymptotically 
large n and R, we adopt it now as part of the definition of the shell model 
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for arbitrary w and R. At this point it may be verified that the necessary 
conditions (12. 6p both hold if 

(c d n)-^ < 1. (2.11) 
Equation (12.101) is the desired relation between w and R. 

2.3 Maximizing the entropy 

Upon substituting (I2.10p in (12. 4p and denoting the result by p n (d; R) we 
obtain 

logp n {d;R) = -\ogn\ + n\ogpv d {2R) d - pv d {2R) d + n\ogZ n {d) (2.12) 
where 

Z n {d) = l-{l-(c d n)-^r} d 

~ d(c d n)-^, (2.13) 

in which the second line represents the leading order behavior as n — > oo. 

Equation fl 2 . 1 2 [) represents the entropy of the arrangement of seeds and 
still contains R as a free parameter. It is again easy to maximize log p n (d; R) 
expression with respect to R. 

Upon varying the right hand side of (I2.12p with respect to R we find that 
it has its maximum for R = R* where 

2R* =(—) . (2.14) 

The corresponding w* follows from substitution of (12.141) in (12.101) . We are 
now ready to obtain our heuristic result as the probability that maximizes 
the configurational entropy, that is, p n (d) = p n (d; i?*). Substitution of (12.141) 
and fl2TT3|) in fl2TT2l) yields 

\ogp n (d) = — logn! + nlogn — n + nlog 4 c rf n ) _3rT + n ^°S z n(d), (2.15) 
in which 

e n (d) = dr x {c d n)^Z n {d) - 1 

^ (-l) fc / d \ . 2k , 

k=l v 7 
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is a polynomial in (qti) rf- 1 without constant term. This may be rearranged 
to yield the large-n expansion 

logp„(d) = -log[(^;)!] + n log A d + n log [l + e n (d)) + ilog^ + o(l), 

(2.17) 

where is an abbreviation for 

A d = d[\{d-l)c d e]~^ . (2.18) 

Hence we may write 

i A n 

Pn(d) = (^^[l + aCtl + ^l, n^oo, (2.19) 

which is our final result for general d. 

Since l+e n (d) has to be elevated to the nth power, this polynomial cannot 
be included with the o(l) terms when d > 3. Dimension d = 3 is a marginal 
case and e n (3) contributes to the constant prefactor of p n (3). Hence p n (3) is 
given by 

A n 

p n (3) = C-f- [1 + o(l)], n^oo, (2.20) 
n! 

in which C = exp(— I/C3) and A3 = l/(3c3e) are numerical constants. 
2.4 Comments 

The arguments of this section are heuristic. Their validity is best assessed 
by a comparison to the two-dimensional case where exact results are avail- 
able. These strongly suggest that the inverse factorial 1/ (jit)!, which is the 
dominant factor in the large-n behavior in (I2.19P and (12.201) . are exact. They 
also lead us to believe that the functional form of p n (d), that is, an expo- 
nential divided by a factorial, can be trusted. The value of the numerical 
constant Cd and hence of Ad remain, however, undetermined, since Cd appears 
in the theory only as a proportionality constant in the order-of-magnitude 
estimate (12.101) . Finally, the decay of p n (d) with growing n is less fast than 
in dimension d = 2, where p n (2) ~ A^/(2n)\ as shown by (II. ip . 

3 Aboav's law in dimension d = 3 
3.1 The plane T of first neighbors 

For large n the first-neighbor seeds are arranged in a nearly spherical shell. 
Its radius, according to (I2.14p and (12. 2p . is equal to 2i?* = \{Qn/ixp) l l z and 
its width, according to (12.141) and (12.101) . is u>* ~ (pn 2 ) -1 / 3 . This width may 
be set to zero for all considerations of this section, which means neglecting the 
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Figure 4: The plane J-, whose actual width w* ~ re~ 2 / 3 may be set to zero, 
contains the seeds that are first-neighbors to the central seed (the latter is located 
at a distance 2R* ~ n 1 / 3 below the plane and is not shown). S\ and S2 are second- 
neighbor seeds. V\ is the paraboloid of all points equidistant from S\ and from T . 
The 'prism' is the Voronoi cell of an arbitrarily selected first-neigbor seed Fq. Each 
vertical face of this prism lies in a plane equidistant from a pair of points in T. The 
intersection of the prism with T is the boundary of the two-dimensional Voronoi 
cell of Fq in T . The prism intersects the parabola V\ according to a polygon which 
is the projection of this boundary. It defines the face that the Voronoi cells of Fq 
and Si have in common. In the downward direction the Voronoi cell of Fq ends 
by a face (not shown) in common with the central cell. 
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Figure 5: The piecewise paraboloidal surface T is the locus of points equidistant 
from the plane of first-neighbor seeds, J 7 , and the set of second-neighbor seeds {Si} 
constitutes a piecewise paraboloidal surface called T. All contributing paraboloids 
are circular and have their axes perpendicular to J- . In particular, V\ and V2 are 
the paraboloids of points equidistant from the seeds S\ and S2, respectively, and 
the plane T . The plane Q bisects the line segment S1S2 perpendicularly. The 
curved segment AB lies on the parabola along which V\, V2, and Q intersect. The 
Voronoi cell of the first-neighbor seed F\ happens to intersect V at the intersection 
AB between two paraboloids. As a result this Voronoi cell has not one, but two 
faces at its upper end. 
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small random deviations of the radial coordinates of the first neighbors. The 
surface area of the sphere being S = 4:ir(2R*) 2 = 7r 1 / 3 (6n/p) 2//3 , the typical 
interseed distance between the first-neigbor seeds is £ ~ (S/n) 1 ^ 2 ~ (p 2 n) -1 / 6 . 
On the scale of the interseed distances we may therefore consider the shell 
as a flat surface that we will denote by T and also refer to as a 'plane'. This 
situation has been represented in figure HI 

In the limit of large n the first-neighbor cells become very elongated 
prism-like objects, as already begins to be apparent in figured] (snapshots of 
realistic two-dimensional many-sided cells with n as high as 1500 are shown 
in reference [15]). Each first-neighbor cell has one face in common with the 
central cell. With the width of T set to zero, the faces between adjacent 
first neighbors are perpendicular to T and define the sides of a prism around 
each first meighbor. These prisms intersect the plane T according to a two- 
dimensional cellular structure. The typical cell area in T is a = S/n ~ C 2 ~ 
(p 2 n) _1//3 . The seeds are not uniformly (Poisson) distributed in T but will 
effectively repel each other; nevertheless, just for topological reasons, the 
cells in the plane T have an average of exactly six neighbors. 

3.2 First and second neighbors 

We consider now the faces between the first and second-neighbor cells. The 
second-neighbor seeds are marked Si, S2, ■ ■ ■ in figures H] and [51 There is 
no restriction on their positions as long as they stay out of the sphere of 
radius and the typical distance p~ 1//3 between them is independent of 
n. We denote by T the surface of points that are equidistant from T and 
from the set {S{} of second neighbors. Since for n — > 00 the first-neighbor 
seeds become infinitely dense in J 7 , in that limit T is a piecewise paraboloidal 
surface. The paraboloids join along lines of intersection ('seams') that are 
segments of parabolas. For example, the curve AB in figure [5] lies on such 
a parabolic seam. For the considerations that follow it will be convenient to 
project T onto T . The set of parabolic seams of T will project onto T as a 
two-dimensional cellular net of trivalent vertices, connected by segments of 
parabolas. We will refer to the cells of this network as 'supercells' in order 
to distinguish them from the 'ordinary' cells (discussed above) due to the 
intersections of the prisms with T . The typical supercell area will be of 
order ~ n° as n — > 00. Since the radius of the spherical shell behaves as 
R* ~ n 1 / 3 , it is well approximated by the flat surface T also at the scale of 
the supercells. 

We analyze now, within the plane JF, the intersection of the net of su- 
percells with the ordinary cells. In the limit n — > 00 the fraction of ordinary 
cells not intersected by a segment of the supercell net will tend to unity. For 
reasons that will become clear just below we denote this fraction of ordinary 
cells by /§. The prism (three-dimensional cell) that encloses a cell of this 
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type, will therefore be bounded below by the central cell and above by a 
single second-neighbor cell. Since (for mere topological reasons) such cells 
are adjacent to, on average, six other first-neighbor cells, their total number 
of neighbors is eight. 

There are, however, two special types of ordinary cells: (a) those in- 
tersected by a perimeter segment of a supercell; and (b) those containing 
the vertex where three such perimeter segments join. In figure |5] the cell 
of seed Fx is an example of a special cell of type (a). These two special 
types of cells will represent fractions of all ordinary cells to be denoted fg 
and /io, respectively. A counting similar to the one above easily shows that 
the corresponding types of three-diemensional cells have 9 and 10 neighbors, 
respectively. 

For n — > oo the fractions fg and /io will vanish, and we will now determine 
exactly how. Since the supercells are of linear dimension ~ n° and the 
ordinary cells of linear dimension supercell will contain ~ n ' 

ordinary cells. Only a finite number of these (on average six) will be located 
on the vertices of the supercell, and therefore we deduce that /i = a 2 n~ 1 / 3 
as n — > oo, where a 2 is a numerical constant. The parabolic segments of 
a supercell perimeter are straight lines at the scale of the ordinary cells. 
Typically such a perimeter segment will therefore intersect ~ n 1 / 6 ordinary 
cells. Hence fg is of order n~ 1 ^ . Assuming an expansion in powers of rT x ^ 
we will write fg = a\rT 1 ^ + bxn~ 1 ^ 3 , where ax and bx are numerical constants. 

To order n -1 / 3 we have that fs + fg + fxo = 1- Hence to this order the 
average number of neighbors m n (3) of a cell with n neighbors is given by 

m„(3) = 8/ 8 + 9/ 9 + 10/ 10 . (3.1) 

Substituting the above expressions for the f n yields 

m n (3) = 8 + k in ~ 1/6 + k 2 n~ 1/3 + . . . , (3.2) 

in which kx = «i and k 2 = bx+2ax- The constants kx > and k 2 are unknown. 
Equation (13.21) is the two-dimensional counterpart of (11. 2p . It shows that in 
dimension d = 3 Aboav's linear law cannot hold for n asymptotically large. 
This law therefore is necessarily an approximation (and possibly a very good 
one) in the experimentally accessible window of n values. Below we will 
study the deviation from Aboav's law numerically. 

4 Comparison to simulation data 

In the preceding sections we derived results that are asymptotic in n, whereas 
data are mostly in a certain range of small n. In earlier work it turned 
out [11], however, that the asymptotic expression for the nearest-neighbor 
correlation m n (2) gives a good approximation to the two-dimensional data for 
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all n values. We may therefore hope that equation (13.21) will similarly provide 
a good fit to the three-dimensional data. Only few such data exist. The ones 
most relevant are due to Kumar et al. [T7]. These authors determined by 
simulation, among several other quantities, the probability p n (S) and the 
correlation m n (3) in the range 10 < n < 22. Their p n (3) data, not shown 
here, peak at n = 15. We have reproduced their m n (3) data in our figure 
[6j Following Earnshaw and Robinson [18], we plot function of 

1/n (rather than m n or nm n as a function of n). In the m n versus 1/n 
plot Aboav's law again corresponds to a straight line, but deviations from 
linearity are easier to detect. 

Within the measurement window m n (3) is clearly seen to decrease with 
n, but only from around 16.4 to 16.1. Kumar et al. fitted this behavior by 

m n (3) = 16.57 -0.02ra, (4.1) 

represented by the dashed line in figure [6J This relation obviously cannot be 
asymptotic. In a later analysis of the same data, Fortes [TS] proposed to fit 
them by Aboav's law, namely 

4.45 , 
m„(3) = 15.95 + — . (4.2) 

It is shown as the straight dotted line in figure EJ 

We wish to compare these two earlier fits to our theoretical functional 
form, equation (13.21) . To that end we choose the constants k% and hi such 
that in n = 15 (where p n is maximum) our curve and the dashed fit produce 
the same values of m n (3) and its n-derivative, that is, rn.xs(3) = 16.27 and 
m' 15 (3) = —0.02. This leads to k\ = 23.15 and k 2 = —15.96. The result is the 
solid curve shown in the figure. We emphasize that this procedure involves 
the additional assumption that it is correct in the finite-n regime to use the 
asymptotic expression (13.21) with all terms beyond order n _1//3 discarded. 

Returning now to figure O we observe that the simulation data scatter too 
much to be able to unambiguously distinguish between the three curves. The 
considerations of this section point to a most interesting question: can one 
establish by Monte Carlo simulation the presence of the downward curvature 
in the m n versus 1/n plot in three dimensions? Curvature, although not 
proving (I3.2p . would at least rule out Aboav's law. Simulations at least 
an order of magnitude larger than the existing ones will be necessary; this 
however is within present-day machine capacity. 

5 Conclusion 

We have considered Poisson-Voronoi diagrams in spatial dimensions d higher 
than two. We obtained analytic expressions for (i) the facedness (or : neigh- 
bor number) probability p n {d) and (ii) the two-cell correlation m ra (3), both 
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Figure 6: Dots: data points by Kumar et al. |17| . Straight dotted line: Aboav's 
law with parameters as in (|4.2|) due to Fortes [T3] . Dashed curve : nonasymptotic 
fit (|4.ip due to Kumar et al. Solid curve: our theoretical equation (|3.2p for the 
values of &i and &2 given in the text. 

valid in the limit of asymptotically large n. We conclude that Aboav's law 
cannot be strictly valid in dimension d — 3, although it may be a very good 
approximation in the regime most easily accessible to experiment and sim- 
ulation. These results rest on heuristic arguments developed in analogy to 
reasoning previously shown [9] to be valid in two dimensions. They cannot 
be considered as mathematically proved, but their value derives from the fact 
that this is the only theoretical work so far in this direction. We believe that 
confirmation of the failure of Aboav's law in three dimensions is within the 
reach of Monte Carlo simulations that are possible today. 
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